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ABSTRACT 


The solution of the Hallen integral equation for a resonant cylindrical an- 
tenna is presented in a simple trigonometric form. The method of solution is 
discussed. A Romberg integration scheme, used in the antenna calculations, is 
given. The computations are discussed and a sample of results is shown. 
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CURRENT AND DRIVING POINT 
IMPEDANCE COMPUTATIONS FOR 
LONG RESONANT ANTENNAS 

I. DESCRIPTION OF THE EQUATION 

Finding a representation for the current on long cylindrical antennas that is 
both simple and in good agreement with more complicated formulations is diffi- 
cult because such a representation is an approximate solution to the Hallen in- 
tegral equation for the current, I 2 (z): 

I* I z (z')K(z, z')dz* = - [c cos(/3 o z) + yV sin(/3 0 |z| )] (1) 

J- h 

where 


K(z.z') = R" 1 exp(-j/3 o R) 


with 


R = ((z - z 1 ) 2 + a 2 ) 2 

£ o - 120 ohms, /3 q - 2rr /\ , 


a - radius of the antenna 
f3 o h - n77/2 n = 1,3, .... 


2h - length of the antenna in wavelengths 


The antenna is center-driven by the delta -function generator, V q . 
This is an integral equation of the first kind which has the form 


r 

a 


f(z' )K(z,z')dz' - h(z) , 


1 



where K and h are known, f unknown, functions of z; as opposed to integral 
equations of the second kind which have the form 



f(z')K(z, z' )dz' +h(z) 


f(z) . 


where again h and K are known, f unknown. 

From the viewpoint of the physicist equations of the latter type are easier 
to solve because a specific formulation for f does not have to be proposed be- 
fore numerical methods of analysis can be meaningfully employed to obtain the 
desired solution. Quite the opposite with equations of the former type: experi- 
ence shows that a model for f , which adequately describes the physical situation 
as known from experiment, must be devised before the parametric values ob- 
tained in the approximate solution will yield a good fit to the experimental data. 
That is, no amount of mathematical analysis, however excellent can make a poor 
model for f yield good results. 

Using experimental data as a guide, King and Sandler (3) have derived a 
simple trigonometric current representation for long resonant antennas. The 
antenna current consists of a shifted sine and cosine distribution with coefficients 
A and B. It has the following final form: 


I 2 (z) “ “ j A cos [/3 o (z-m o )] + B cos (/3 o z ) + g(z, m o ) (2) 

where 

A = 2tT / O ^CR ( m o ) sin (A, m o )] 


B = i A [cp ci (h) 


+ <Pc ‘ ^CR 


(m o ) cos {/3 0 (h-m o )} /qp C2 (h) 


(4) 


g(z,m o ) - ~ jA cos {/3 0 (h-2m o )} cos(/3 o z)/cos {/3 0 (h-m o )} 


(5) 


- cos {/3 o (z-m o )} J h-m o <z<h 


g(z,m o ) - 0 


0 < z < h-m 

— o 


2 


^CR ( m o ) " REAL £ cos (/3 o m o ) C a (h,m o ) - sin (/3 o m o ) S a (h, m Q )J 

(6) 

q>d(h) - cos (/3 o 

m o> C a( h * m o- h ) + Sin ^o m o> S a( h " m o' h ) 

(V) 


cp C2 (h) = C a (h,h) 

(8) 


= C a (h,h) - C a (h-m o ,h) 

(9) 

c a (h, z) 

= f cos (yS o z')K(z, z 1 ) dz‘ 

X h 

(10) 

S a (h, z) 

- f sin (/3 o | z' | K(z, z 1 ) dz‘ 
J-h 

(11) 


The phase shift m o must be determined from a transcendental equation given 
by 


Cos (/3 0 m 0 ) - sin (£ 0 m o ) - 


- REAL 



( 12 ) 


(h) + q?c " <P C R( m o) cos P o ( h 


_m o>} 


Once the value of m is calculated, A and B can then be determined from 

O 7 

Equations (3) and (4). 

The first model for the current gave a poor fit to the experimental data be- 
cause (5) was not well formulated. Consequently, this equation was modified so 
that the final representation for I z (z) compared excellently with both experi- 
mentally obtained and theoretically calculated currents/ 3 ) 
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n. METHOD OF SOLUTION 


The function f (m ) is defined by 

f (m) = cos (/3 o m) - sin (/? 0 m) 

(13) 

+ REAlT q?c2 ( h ){ <Pci( h ) + “ ^*cr cos { ( h “ m )} • 


Clearly, a zero of (13) is a solution of (12). From physical considerations m o is 
such that 0 < m < .25. 

O 

The Newton-Raphson method with suitable restrictions to insure that nr was 
in the first quadrant for i = 1, 2, N, was chosen to find m o . This is an itera- 
tive procedure of the form 

m i + 1 = nii — f (m i ) / f ' (m i ) , i = 1,2, . . . , N . (14) 


The recursion formula is given by setting <£(m) = m - f(m)/ f'(m). The resulting 
relation m i + 1 - 0(nr) is the recursion formula. The sequence {m n ) converges 
to the fixed point, m Q =0 (m o ), of <£ provided that m x is sufficiently close to m o , 
that $ (m) is continuous in an interval around m o , and that \cf> (m o ) | < l.CS) In 
the Newton-Raphson method 4> (m) = f (m) f" (m)/( f'(m)) 2 , so that <f> (m o ) = 0 
provided that f is twice continuously differentiable and that f'(m o ) ^ 0. 

The analytic expression for f'(m) would be extremely difficult (and tedious) 
to find and equally cumbersome to programe for evaluation. Thus the derivative 
of f (m) was approximated by the difference formula 


f' (m) 


f (m + Am) - f (m) 
Am 


(15) 


where Am was sufficiently small. 

All integrals were evaluated numerically. A set of recursion formulas, 
called a Romberg scheme for Simpson sums, was used. The method is based on 
repeated interval halving until the desired degree of accuracy is reached. 
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m. ROMBERG SCHEME FOR NUMERICAL INTEGRA TION ( 5 ) 


This section outlines a quadrature scheme, based on interval halving, which 

evaluates numerically the definite integral J* f (x) dx, where f (x) is a continuous 

function that can be calculated by some means at any point x. (To facilitate dis - 
cussion the interval [a,b] is normalized to [0,1] by a suitable transformation 
on x). First a recursion formula for trapezoidal sums is described. The 

trapezoidal rule for mesh size, h, gives the approximation, T, to f (x)dx: 

Jo 


T 



2 {f (h) + f (2h) + ... + f ( 1 - h) } + f (1)" 


( 4 ) 


The recursion expresses the i th trapezoidal sum, T., of mesh-size, h = (■§) l * 
as j the sum of T. x and the mean of the values f (x) when x takes on the values 
(2j - l)/2 iml , j = 1, 2, ...,2 i ‘ 2 . Then a second recursion formula is derived 
to express the i th Simpson sum, S., of mesh-size, h = (y) 1 , as a linear com- 
bination of T i + 1 and Simpson's rule approximates f (x)dx by the sum, S, 

Jo 

S = | [f (o) + 4 { f (h) + f(3h) + . . . + f(l-h)} 


+ 2 {f (2h) + ... + f (1 - 2h)} + f (1)] , (4) 


for the mesh, h. 

This recursive method has the outstanding advantage that as one decreases 
the mesh size from (-|) n to (j) n + 1 , one need compute only 2 n new points, instead 
of the 2 n+1 + 1 points usually computed. Thus, computation time is saved. 
Furthermore, the mesh-size can be made a function of the convergence criterion. 

In Figure 1, the area of trapezoid ACDF, denoted by T (|), being the mid- 
point of its base along the X-axis , is a rough approximation to the area under the 
curve y = f(x) in the interval [0,1]. Define 

Ti = T(|) = j[f(o) + f(l>] (16) 
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Figure 1. 


0 A 'A 

Hopefully the sum of the area of trapezoid ABEF , denoted by T (-£•); and of BCDE, 
designated T(|), is a better approximation, T 2 , to the desired area. Thus, if 
T(i), T(|-) and T 2 are defined by the following relations 

T(i) = \ [f (o) + f(i)] 

T( | ) = + f ( x >] (17) 

T 2 = T(i) + T(|) = H T i + f <2>] 

it is hoped that 

| f f (x) dx - T 2 

Jo 

In the same manner the area under y = f (x) may be approximated by the four 
trapezoidal areas: 

T(f) = I D(0) + f(?)] 


4 

^0 


f (x) dx 


6 





T (| ) = i[f(i)+ f(i)] 


( 18 ) 


T(|) = i [f(4) + f(i>] 

T(|) n t [f(D + f Cl)] 


The sum of these areas in (18) is T, where 




A T + 
2 A 2 ^ 


f(-b + f(i) 


Using the first inductive principle one can prove that if 


T x = | [f CO) + f Cl)] 


(19) 


then 


T 



i - 2, 3, . . . . 


where 


( 20 ) 


N = 2 i_2 , M = 2 i_1 


For the unnormalized interval [a,b] 

T x = [f (a) + f(b)] 


and 


T. = 


b~ a 


T i-l 
b- a 


j f(at {^} {b . a}) 

i=i 


/N 


i = 2,3. 


(20’) 
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where 


N = 2 i ~ 2 , M = 2 


i - 1 


It is obvious that 


1 i m ’T 


r 

•'a 


f (x) dx 


( 21 ) 


The following paragraphs will now establish a linear relationship between 
elements of the sequence ^ °f successive approximations to the (definite 

integral I f (x) dx using Simpson's rule and the elements of the sequence {T. ]• 

of successive approximations using the trapezoidal rule. 


In mechanical quadrature f (x) is replaced by an interpolating formula inte- 
grated between the limits, a and b, (in this discussion 0 and 1), df the definite 
integral. For Simpson's rule the interpolating formula is a polynomial of de- 
gree 2, denoted P 2 (x). Therefore, a linear combination, S 1 , of T x and T 2 , is 

T 1 

required which yields exactly the value of the definite integrals, I P 2 (x)dx. 
For example, if f (x) = x 2 , 


I’ 

** r\ 


x 2 dx = j -- 


pTj + qT 2 


( 22 ) 


with p + q = 1. Since Tj = T 2 = |, the following system of linear equations 
is solved: 


p t q 


1 


p/2 + 3 


s'! ~ 


(23) 


obtaining p 


i _ 

3 * 


= so that 


S z = (4T 2 - Tj ) / 3 


(24) 
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Geometrically, Simpson's rule replaces f (x) by parabolas (the trapezoidal 
rule, by chords), as follows: for a mesh of j, a parabola, denoted by P 2 . ^ ( x) * , 
is passed through the points (0, f (0) ) , ( i, f (-1) ) , (1 , f (1) ) , which replaces f ( x ) in 

the interval [0,1]. Thus, I f (x)dx is approximated by I P 2 ,^(x)dx. If the 

■'o •'o 


interval [0,1] is divided into quarters, f(x) can be approximated by the sum of two 
parabolas, P 2 , /4 (x) through the points (0, f(0) ), (}, f(l) ), (1, f(|) ); and P 2 (x) 
through (7, f(y) ), (f> f ( -§-) ), (1, f(l) ). Integrating this sum one obtains, using (24) 

and (20’) for intervals of length \ , the approximation, S 2 , to f f (x) dx, 


(4T 3 -T 2 )/3 


(25) 


since 


Ti = |[f(o) + f(i)] 


T i = i[ f (i) + f (!)] 


V 2 = i[2Ti + f(i)] 


T-J = i[2T" + f(|)] 


s; = (4 T'-t;>/3 


s" = (4T 2 - r; )/3 


s 2 = s; + sj = [t; + t; + f(|) + f(|>] /3 = (4 t 3 - t 2 )/3 


where 


' applies to sums on the interval [0,1] , 

" to sums on [ y, 1] . 

By dividing the interval [0, 1] into eighths, it can be shown that 


S 3 = (4T 4 -T 3 )/3. 




2 ,m 


(x) is a parabola which replaces f(x) in the sub-interval with midpoint, m. 
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Generally, the recursion relation is of the form 


S. 


(4T l + 1 - T i)/ 3 - 


( 26 ) 


The sequence {S n } is a sequence of successive approximations to 
using Simpson's rule with mesh-size h = h n = 2~ n . 

To utilize the Romberg method in evaluating integrals of the form 



f (x) dx, 


f (x) dx, 


a S/360 Fortran IV program was written. It is available from the Goddard Pro- 
gram Library. ^ 2 ^ This program is especially recommended in cases where an 
optimal, but unknown, mesh-size is needed in order to approximate an integral 
to a pre -determined number of significant figures. Test cases indicate that up 
to six figures of accuracy can be obtained using this program. 


IV. DISCUSSION OF COMPUTATIONS AND RESULTS 

The computations were done on the IBM 7094-7040 DCS and the IBM S/360 
ASP systems at Goddard Space Flight Center. The IBM S/360 Fortran IV, 
Version H, program of the King -Sandler theory is documented and available from 
the Goddard Program Library. ( 1 ) 

To evaluate the functions containing integrals, (10) and (11) were split into 
their real and imaginary parts. Each of the resulting integrals was then approxi- 
mated separately using Program No. M00077 in the Goddard Program Library. 
This program utilizes the method discussed in Section III. 

Successive approximations, S. and S i + 1 , to an integral, I, were tested for 
convergence as follows: 

given e , if 


(Si - S i*,)/Si I < £ 


(27) 


then 


I 


S 


i + 1 


For all of the results given in this paper e = 10" 3 was chosen because: 
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1) if e > 10 -3 , the m o obtained in the solution of (13) differed significantly 
from the m o for e = 1 0“ 3 ; 

2) the m o obtained for e < 10“ 3 agreed to three figures with that for 
e = 10" 3 , but the computational times differed greatly. 

In the Newton-Raphson iterative procedure zero was approximated by 10 -3 ; 
that is, if f (m i ) < 10~ 3 , m o = m i the desired zero of (13). Once again this cri- 
terion was established empirically: such a tolerance seemed to give an optimal 
trade-off between time and accuracy. Similar considerations indicated that 
Am = 10" 3 be used in (15). 

Each iteration required two evaluations of (13) which, in turn, was a function 
of the integrals given by (6), (7), (8), (9). On the IBM S/360 ASP, 57.5 seconds 
were required to evaluate (13). An average of 5.3 iterations was needed for the 
convergence of the sequence {nr } to m o . 

However, a good estimate for A shortened the number of iterations to 
one or two; while a bad estimate caused an end of case condition to occur. To 
facilitate the selection of m t A for a given h A, a plot of m o A as a function of 
hA was made for completed cases. Figure 2 shows such a graph for Q. - 20. 
Since Line I appeared to be an upper bound for the points (h A, m o A); and Line II 
a lower bound, (hA, m x A) less than (or equal to) the corresponding point on 
Line I and greater than (or equal to) that on Line II was chosen as the starting 
point. 

Table I summarizes the King-Sandler results for fl = 20. For a given half- 
length, h A, the phase shift, m o A, the driving point admittance, Y o , the driving- 
point impedance, Z o , and the coefficients, A and B, calculated using (3) and (4), 
respectively, are tabulated. All computations are to four significant figures. 

The parametric values of h A, m o A, A and B can then be used in (2) to compute 
the current at any point along the antenna. 

The good agreement between the King-Sandler approximate model and the 
experiments of Altshuler and theories of Wu and King-Middleton for selected 
values of h and fi indicates that this model is an adequate description of cur- 
rents on long resonant antennas. 
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Table I 


Resonant Antenna Parameters 
A, B, Y 0 in milliohms, Z c in ohms 


n = 20 


hA 

m <A 

Yo 

z 0 

A 

B 

0.25 

0.0455 

10.08 - j5.23 

78.2 + j40.6 

3.19 

10.08 - j5.23 

0.75 

0.0719 

7.70 - j3.04 

112.3 + j44.4 

2.23 

7.70 - jl.03 

1.25 

0.0853 

6.86 - j2.42 

129.7 + j45.8 

2.00 

6.86 - j 0.70 

1.75 

0.0941 

6.37 - j2.09 

141.7 +j46.6 

1.90 

6.37 - j0.52 

2.25 

0.1005 

6.03 - jl.88 

151.1 +j47.2 

1.84 

6.03 - j0.40 

2.75 

0.1056 

5.78 - jl.73 

158.9 + j47.6 

1.80 

5.78 - j 0.31 

3.25 

0.1098 

5.57 - jl.62 

165.6 + j48.0 

1.78 

5.57 - j0.24 

3.75 

0.0521 

8.68 - j5.29 

84.0 + j51.2 

3.36 

8.68 - j 2 . 11 

4.25 

0.0526 

8.62 - j5.27 

84.4 + j51.6 

3.36 

8.62 - j2.09 

4.75 

0.1189 

5.14 - jl.38 

181.5 + J48.1 

1.76 

5.14 - j0.10 

5.25 

0.1212 

5.03 - jl.33 

185.9 + j49.1 

1.75 

5.03 - j0.06 

5.75 

0.1233 

4.93 - jl.28 

190.0 + j49.4 

1.75 

4.93 - j 0.03 

6.25 

0.1251 

4.84 - jl.24 

193.9 + j49.6 

1.75 

4.84 + jO.OO 

6.75 

0.1267 

4.76 - jl.20 

197.5 + j49.7 

1.76 

4.76 + j0.03 

7.25 

0.0813 

6.90 - j3.02 

121.6 + j53.2 

2.40 

6.90 - j0.92 

7.75 

0.0553 

8.36 - j5.18 

86.5 + j53.6 

3.36 

8.36 - j2.02 

8.25 

0.0555 

8.33 - j5.17 

86.7 + j53.8 

3.36 

8.33 - j2.01 

8.75 

0.0823 

6.83 - j3.01 

122.7 + j54.1 

2.42 

6.83 - j 0.91 

9.25 

0.1333 

4.45 - jl.05 

213.0 + j50.5 

1.77 

4.45 + j 0.13 

9.75 

0.1344 

4.39 - jl.03 

215.7 + j50.6 

1.78 

4.39 + j0.15 

10.25 

0.1354 

4.35 - jl.01 

218.3 + j50.7 

1.78 

4.35 + j0.17 

10.75 

0.1364 

4.30 - j0.99 

220.9 + j50.8 

1.79 

4.30 + j 0.18 

11.25 

0.1373 

4.26 - j0.97 

223.3 + j50.9 

1.79 

4.26 + j0.19 

11.75 

0.0573 

8.17 - j5.10 

88.1 + j55.2 

3.36 

8.17 - j 1.96 

12.25 

0.0576 

8.15 - j5.09 

88 .2 + j 55.2 

3.36 

8.15 - jl.95 

12.75 

0.1396 

4.14 - j0.92 

230.2 + j51.2 

1.81 

4.14 - j 0.23 

13.25 

0.1403 

4.10 - j0.91 

232.4 + j51.3 

1.81 

4.10 - j 0.25 

13.75 

0.1410 

4.07 - j0.89 

234.5 + j51.4 

1.82 

4.07 - j 0.26 

14.25 

0.1417 

4.03 - j0.88 

236.6 + j51.5 

1.82 

4.03 - j0.27 

14.75 

0.0976 

6.00 - j2.38 

144.1 +j57.2 

2.25 

6.00 - j 0.54 

15.25 

0.0857 

6.59 - j2.97 

126.2 +j56.9 

2.47 

6.59 - j0.85 

15.75 

0.0592 

7.69 - j5.01 

89.9 + j56.5 

3.34 

7.97 - jl.90 
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